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Abstract. A phenomenological systems approach for iden- 
tifying potential precursors in multiple signals of different 
types for the same local seismically active region is proposed 
based on the assumption that a large earthquake may be pre- 
ceded by a system reconfiguration (preparation) at different 
time and space scales. A nonstationarity factor introduced 
within the framework of flicker-noise spectroscopy, a statis- 
tical physics approach to the analysis of time series, is used 
as the dimensionless criterion for detecting qualitative (pre- 
cursory) changes within relatively short time intervals in ar- 
bitrary signals. Nonstationarity factors for chlorine-ion con- 
centration variations in the underground water of two bore- 
holes on the Kamchatka peninsula and geacoustic emissions 
in a deep borehole within the same seismic zone are studied 
together in the time frame around a large earthquake on Oc- 
tober 8, 2001. It is shown that nonstationarity factor spikes 
(potential precursors) take place in the interval from 70 to 
50 days before the earthquake for the hydrogeochemical data 
and at 29 and 6 days in advance for the geoacoustic data. 



1 Introduction 

Earthquake prediction in the time frame of several months 
to less than an hour before the catastrophic event, which is 
often referred in literature as "short-term" prediction, has 
been a subject of extensive research studies and controver- 
sial debates both in academia and mass media in the past 
two d ecades (|Geller||1997 



1997 



Wyss et al. 



Geller et al. 

[T997||Uyeda et al.| |2009[ [Cicerone et aLl ' |2009l l. One of the 
key areas in this field is the study of earthquake precursors, 
physical phenomena that reportedly precede at least some 
earthquakes. The precursory signals are usually grouped into 
electromagnetic, hydrological/hydrochemical, gasgeochemi- 
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cal, geodetic, and seismi c ([Geller 1997 [Hartmann and Levy 



2005 1 |Uyeda"etaLl [20091 [Cice rone et al.' '20091 Electromag 
netic precursory signals are further classified into the signals 
believed to be emitted from within focal zones, such as tel- 
luric and magnetic field anomalies, and radio waves over epi- 



central regions ( Uyeda et al. 2009). The localized changes in 



electric and magnetic fields that reportedly accompany some 
seismic events span a wide range of frequencies, including 
ULF, VLF, ELF and RF fields, and were observed in the time 
frame from 2-3 years to dozens of minutes prior to an earth- 
quake ( [Cicerone et"aLl [2009] [Uyeda et a\\ [2009^. Hydro- 
logical/hydrochemical precursory signals include water level 
or quality changes in the weeks, days, or hours prior to a 
number of earthquakes, groundwater temperature changes, 
and variations in the concentrations of dissolved ions like 
chlorine or magnesium usually in the time frame of months 
to days before an earthquake (Hartmann and Levy 2005 
Cicerone et al. 2009; Du et al. ,2010| [ Gasgeochemical 



precursory signals comprise numerous anomalous gas emis- 
sion observations, the majority of which were reported for 



the concentration of radon gas in the earth (Hartmann and 
Levy] [2005) [Cicerone et~aL| [2009|. More than 100 stud 



ies show that changes in radon exhalation from the earth's 
crust precede a number of earthquakes by months, weeks, 
or days ( [Cicerone et al.| [2009[ ). Geodetic signals mostly in- 
clude surface deformations (tilts, strains, strain rate changes) 
over distances of tens of kilometers that precede some ma- 



jor earthquakes by months to days ( Cicerone et al. 2009 1. 



Seismic precursory signals encompass foreshocks that typi- 
cally take place less than 30 days before the main shock and 
high-frequency (acoustic emission) and very low-frequency 
precursory signals that are not detected by conventional seis- 



mographs (Dimle and Jordan, 1994, Reasenberg 1999 Gor- 



dienko et al. 2008 Gavrilov et al. 2008 1. Another promis- 
ing type of possible precursory signals is anomalous animal 
behavior for very short time frames (within 2-3 days, usu- 



ally hours) prior to a large seismic event ( Kirschvink 2000 
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lYokoi et al.||2003i|Li et al.[|2009l ). 

Despite the large number of earthquake precursors re- 



ported in literature, most of which are summarized by 'Hart- 



|mann and Levy] ( |2005[ ); [Cicerone et al.| (2009 ), an Interna 
tional Commission on Earthquake Forecasting for Civil Pro- 
tection concluded on October 2, 2009, "the search for pre- 
cursors that are diagnostic of an impending earthquake has 
not yet produced a successful short-term prediction scheme" 
( |ICEFCP| |2009}. The reports of the International Associa- 
tion of Seismology and Physics of the Earth's Interior con- 
tain similar findings (Wyss and Booth 1997| l. The lack of 
confidence can be attributed to several reasons. First, some 
fundamental aspects of many non-seismic signals, for exam- 
ple, lithosphere-atmosphere-ionosphere coupling and propa- 
gation of high-frequency electromagnetic signals in the con- 
ductive earth, are unresolved, and many of the proposed 
physical models are questionable (Uyeda et al. 2009[ l. Sec- 
ond, the experimental data on precursory signals are often 
limited to few earthquakes and few measurement sites, they 
frequently contain gaps and different types of noise (Hart-| 



mann and Levy 2005 Cicerone et al. 2009 Uyeda et al 



2009 1. Third, different techniques of identifying the anoma- 



lies are used for different signals or even in different studies 
for the same signal. In some cases, the anomalous changes 



are determined by analyzing the signals themselves (Hart 



mann a nd Levy| |2005t |Uyeda et~ar. 2009 ; Cicerone et al. 



|2009jl, while in other cases they are identified by studying 
the derived statistics or functions, such as Fisher information 
or scaling parameters ( Telesca et al.| 2009a|b| l. Moreover, 
seasonal changes and instrumentation or other background 
noise often need to be filtered out prior to the identification 
of precursors. 

In view of the above three problems, we believe that earth- 
quake precursor research can be advanced by employing a 
phenomenological systems approach to the analysis of sig- 
nals of different types in the same local geographic region. 
We assume that a large earthquake may be preceded by a sys- 
tem reconfiguration (preparation) at different time and space 
scales, which manifests itself in qualitative changes of vari- 
ous signals within relatively short time intervals. For exam- 
ple, such anomalous hydrogeochemical signals may be ob- 
served months to weeks before the impending earthquake, 
anomalous geoacoustic emissions - only days prior to the 
event, and anomalous behavior of animals - only hours be- 
fore the catastrophe. In order to test this approach and iden- 
tify different signals that may be related to a specific large 
seismic event, one needs to have a standard criterion or a 
set of standard criteria to detect signal anomalies in virtu- 
ally arbitrary signals. In this study, we will use a nonsta- 
tionarity factor introduced within the framework of flicker- 
noise spectroscopy (FNS), a statistical physics approach to 
the analysis of time series ( Timashev and Polyakovj |2007[ 
[Timashevl |2007| [Timashev et al. 2010b| l. This dimen- 
sionless criterion is practically independent from the indi- 
vidual features of source signals and is designed to detect 



abrupt structural changes in the system generating the sig- 
nal, which makes it a promising candidate to be one of the 
standard criteria. The nonstationarity factor was previously 
used to detect precursors in electrochemical and telluric sig- 
nals recorded in the Garm area, Tajikistan prior to the large 



1984 Dzhirgatal earthquake (Descherevsky et al. 2003 Vs 



tovsky et al. 2005 i, geoelectrical signals at station Giuliano, 



Italy prior to several 2002 earthquakes ( Telesca et al. 2004 1 



and ULF geomagnetic data at Guam prior to the large 1993 
Guam earthquake ( Hayakawa and Timashev[[2006[|Ida et al.j 
j2007| l. Other approaches to identifying precursory features in 
earthquake- and volcano-related signals, which are based on 
different monlinear analysis techniques, were discussed by 
Telesca et al. (2010] |2009a|b| ; [Telesca and LovaUo] ( |2009l l; 
[Telesca et al.| ( |2008l . 

In this study, we consider a combined analysis of two 
different types of signals, hydrogeochemical (sampling fre- 
quency of 3 to 6 day^^) and geoacoustic (sampling fre- 
quency of 1 min^^), recorded on the Kamchatka peninsula, 
Russia. 



2 Nonstationarity factor 

Here, we will only deal with the basic FNS relations needed 
to understand the nonstationarity factor. The approach is de- 
scribed in detail elsewhere ( [Timashev 2006 Timashev and[ 
[Polyakovl [20071 |Timashev[ |2007| [Timashev et aL] [2010b[ ). 
The FNS procedures for analyzing original signal V{t), 
where t is time, are based on the extraction of information 
contained in autocorrelation function 



^iT)^{Vit)V{t + T)), 



(1) 



where r is the time lag parameter. The angular brackets in 
relation ([T]) stand for the averaging over time interval T: 



1 f^/^ 
((...)> = -/ {...)dt. 

^ J-T/2 



(2) 



To extract the information contained in ■iP{t), the follow- 
ing transforms, or "projections", of this function are ana- 
lyzed: cosine transforms (power spectrum estimates) S{f), 
where / is the frequency. 



S{f)^ f ' {V{t)V{t + T))cosi2nfh)dh 

J-T/2 



(3) 



and its difference moments (Kolmogorov transient structural 
functions) of the second order $(^)(r) 



^(^Hr) = {[Vit)-V{t + T)f). 



(4) 



To analyze the effects of nonstationarity in real processes, 
we study the dynamics of changes in (f*^^^ (t) for consecutive 
"window" intervals [tk,tk + T], where k = 0, 1, 2, 3, and 
tk = kAT, that are shifted within the total time interval Ttot 
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of experimental time series (tk+T < Ttot)- The time inter- 
vals T and AT are chosen based on the physical understand- 
ing of the problem in view of the suggested characteristic 
time of the process, which is the most important parameter 
of system evolution. The phenomenon of "precursor" occur- 
rence is assumed to be related to abrupt changes in functions 
<i>'^^)(T) when the upper bound of the interval [tk,tk+T] ap- 
proaches the time moment tc of a catastrophic event accom- 
panied by total system reconfiguration on all space scales. 

The analysis of experimental stochastic series often re- 
quires the original data to be separated into a smoothed and 
fluctuation components. In this study, we apply the "re- 



laxation" procedure proposed by Timashev and Vstovskii 



( 2003 I based on the analogy with a finite-difference solu- 
tion of the diffusion equation, which allows one to split the 
original signal into low-frequency Vii{t) and high-frequency 
Vpit) components. The iterative procedure finding the new 
values of the signal at every relaxation step using its val- 
ues for the previous step allows one to determine the low- 
frequency component VB{t). The high-frequency compo- 
nent Vpit) is obtained by subtracting V«(t) from the original 
signal. This smoothing algorithm progressively reduces the 
local gradients of the "concentration" variables, causing the 
points in every triplet to come closer to each other. Such 
splitting of the original signal V{t) into Vii{t) and Vpit) 
makes it possible to evaluate the nonstationarity factor for 
each of the three functions Vj{t) (J = R, F, or G), where 
index G corresponds to the original signal. 

The FNS nonstationarity factor Cj(tk) is defined as 



Qi-pi ^, T 



AT' 



[Vj{t)-Vj(t + T)fdtdT, 



aT tk+T- AT 



[Vj{t)-Vj{t + T)fdtdT. 



(5) 



(6) 



(7) 



Here, J indicates which function Vj{t) (J = R, F or G) is 
used. Expression s (|6l|7l) are given in discrete form elsewhere 



(Timashev et al. 



|2010a^ . Note that functions $j ^(r) can be 



reliably evaluated only on the r interval of [0, aT], which is 
less than half of the averaging interval T; i.e., a < 0.5. 



3 Experimental data for the Kamchatka peninsula 

The data were recorded in the south-eastern part of the Kam- 
chatka peninsula located at the Russian Far East. The eastern 
part of the peninsula is one of the most seismically active re- 
gions in the world. The area of highest seismicity localized 
in the depth range between and 40 km represents a narrow 
stripe with the length of approximately 200 km along the east 



coast of Kamchatka, which is bounded by a deep-sea trench 
on the east ( |Fedotov et al.||1985| l. 

Specialized measurements of underground water charac- 
teristics were started in 1977 to find and study hydrogeo- 
chemical precursors of Kamchatka earthquakes. Currently, 
the observation network includes four stations in the vicinity 
of Petropavlosk-Kamchatsky (Fig. 1). The Pinachevo sta- 
tion includes five water reservoirs: four warm springs and 
one borehole GK-1 with the depth of 1,261 m. The Mo- 
roznaya station has a single borehole No. 1 with the depth 
of 600 m. The Khlebozavod station also includes a single 
borehole G-1 with the depth of 2,540 m, which is located in 
Petropavlosk-Kamchatsky. The Verkhnyaya Paratunka sta- 
tion comprises four boreholes (GK-5, GK-44, GK-15, and 
GK-17) with depths in the range from 650 to 1208 m. 

The system of hydrogeochemical observations includes 
the measurement of atmospheric pressure and air temper- 
ature, measurement of water discharge and temperature of 
boreholes and springs, collection of water and gas samples 
for their further analysis in laboratory environment. For wa- 
ter samples, the following parameters are determined: pH; 
ion concentrations of chlorine (Cl^), bicarbonate (HCO3"), 
sulfate (S04^), sodium (Na+), potassium (K+), calcium 
(Ca^+), and magnesium (Mg^+); concentrations of boric 
(H3BO3) and silicone (H4Si04) acids. For the samples of 
gases dissolved in water, the following concentrations are 
determined: methane (CH4), nitrogen (N2), oxygen (O2), 
carbon dioxide (CO2), helium (He), hydrogen (H2), hy- 
drocarbon gases: ethane (C2H6), ethylene (C2H4), propane 
(CsHg), propylene (CsHg), butane (C4Hion), and isobutane 
(C4Hioi). The data are recorded at nonuniform sampling 
intervals with one dominant sampling frequency. For the 
Pinachevo, Moroznaya, and Khlebozavod stations, this av- 
erage sampling frequency is one measurement per 3 days; 
for the Verkhnyaya Paratunka station, one measurement per 
6 days. Multiple studies of the hydrogeochemical data and 
corresponding seismic activity for the Kamchatka peninsula 
reported anomalous changes in the chemical and/or gas com- 
position of underground waters prior to several large earth- 
quakes in the time frame from 1987 to 2001 ([Kopylova et al. 



T994i Bella et al. , 1998 ; Bia gi etal.| [2O0O 2006; Khatkevich 
land Ryabinin, |2006) . In this study, we analyze the varia- 
tions of chlorine-ion concentration determined by a titrimet- 
ric method (relative error from 1 to 10%). 

Geoacoustic emissions in the frequency range from 25 to 
1,400 Hz (0.7 level) have also been recorded in the deep G- 
1 borehole of the Khlebozavod station under the supervision 
of V. A. Gavrilov since August, 2000. The data analyzed 
in this paper were obtained by a geophone with crystal fer- 
romagnetic sensors ( jBelyakovj |2000| l. The output signal of 
such a sensor is proportional to the third derivative of ground 
displacement, and the gain slope is 60 dB per decade of fre- 
quency change. The geophone was set up at the depth of 
1 ,035 m, which is enough to reduce anthropogenic noise lev- 



els by more than two orders of magnitude (Gavrilov et al. 
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2008 1. The geophone body was fixed inside the borehole 
casing by a spring. The vertical channel sensitivity of the 
geophone is 0.15 V x s'^/m. The sensitivity of horizontal 
channels is 0.60 V x s^/m. The sensor output signals are sep- 
arated by third-octave band pass filters into four frequency 
bands with central frequencies 30, 160, 560, and 1,200 Hz, 
which is followed by real-time hardware/software signal pro- 
cessing. The value of postprocessed output signal for each 
channel is proportional to the average value of input signal 
for one-minute intervals. More detailed description of geoa- 
coustic emission observations and experimental setup for the 



G-1 borehole is presented elsewhere (Gavrilov et al. 2008 1. 



4 Results 

To illustrate the nonstationarity factor and proposed phe- 
nomenological method, we analyze the hydrogeochemical 
data for chlorine-ion concentrations at GK- 1 (Pinachevo sta- 
tion) and GK-44 (Verkhnyaya Paratunka station) and geoa- 
coustic emissions at the output of geophone vertical fre- 
quency channel with the central frequency of 160 Hz (Z160) 
in G-1 (Khlebozavod station). Chlorine-ion concentration 
time series for GK-1 (Cl-GKl) was selected because it is 
characterized by a unidirectional long-period trend without 
seasonal variations (Fig. 2) and was akeady treated as a pre- 
cursory signal due to a gradual chlorine-ion concentration 
decline down to a local minimum 30 to 60 days before sev- 



eral earthquakes (Khatkevich and Ryabinin 2006) . On the 
other hand, chlorine-ion concentration at GK-44 (C1-GK44) 
is not considered as a precursory signal because it is domi- 
nated by seasonal concentration changes on the background 
of a slowly varying local mean, the minimum value of which 
is reached shortly after the strong earthquake on December 
5, 1997 (Ml =7.0). The Z160 signal was selected from the 
whole set of geoacoustic time series because it contains the 
lowest level of noise (highest signal-to-noise ratio). 

To keep the statistical structure of source time series prac- 
tically intact, the signals were subjected only to minimal 
preprocessing, which included the removal of single-point 
spikes, reduction of the hydrogeochemical time series to uni- 
form sampling intervals using linear interpolations, and ex- 
traction of every 30*^^ point in the geoacoustic time series to 
form a new time series with the frequency of 30 min^^. Then 
the time series Vg (t) were separated out into low-frequency 
Vu{t) and high-frequency Vpit) components, which were 
used to calculate the nonstationarity factors. In evaluating 
Cj [J = R, F or G) for the hydrogeochemical series, aver- 
aging time intervals T in the range from 50 to 900 days were 
used. For the geoacoustic time series, the interval T was var- 
ied from 3 to 20 days. Our analysis showed that the values 
of T equal to 600 and 20 days are most adequate for locating 
precursors in the hydrogeochemical and geoacoustic series, 
respectively. 



Figures 2 and 3 show the variations of Cj for Cl-GKl and 
G1-CK44 together with largest seismic events. It can be seen 
that spikes in Cj precede several large earthquakes. It should 
be noted that the low-frequency component Cr shows the 
most number of precursors for Cl-GKl and high-frequency 
component is most informative for C1-CK44. The first fact 



is in agreement with the study of Khatkevich and Ryabinin 
(2006). The second fact implies that the use of the high- 
frequency component eliminated seasonal changes from the 
analysis and made C1-GK44 a precursory signal. Therefore, 
the FNS nonstationarity factor together with the procedure 
for separating out high-frequency and low-frequency signal 
components can be used to analyze different signals despite 
major differences in their specific features. 

Figure 4 shows a combined analysis of hydrogeochemical 
and geoacoustic variations in the time frame around the Oc- 
tober 8, 2001 earthquake (Mi =6.3, = 24 km, D = 1 34 km 
from Petropavlovsk-Kamchatsky), which was the strongest 
earthquake (based on local magnitude and distance to the 
epicenter) recorded for the whole time interval of geoacous- 
tic observations in the G-1 borehole. Nonstationarity fac- 
tors Cr for Cl-GKl and Cp for C1-GK44 show spikes with 
highest values (precursors) in the time frame from 50 to 70 
days before the earthquake. Cq for G-1 (the signal is a high- 
frequency one by its nature) shows precursors 29 and 6 days 
before the event, which is in agreement with the results re- 
ported by Gavrilov et al. ( 2008| l. In other words, anomalous 
changes in the geoacoustic signal happen closer to the earth- 
quake than in the hydrogeochemical ones, which implies that 
precursory signals of different nature may take place at dif- 
ferent timescales before a large earthquake. 



5 Conclusions 

The above example shows that precursory signals of differ- 
ent types may be observed in the same local seismically ac- 
tive zone at different times prior to a large earthquake, which 
may be attributed to some system preparation preceding the 
seismic event. In the studied case, the qualitative changes 
may be related to a system-wide structural medium reconfig- 
uration at the preparatory phase of the earthquake. 

This study also shows that the FNS nonstationarity fac- 
tor can be used as the standard criterion to detect qualita- 
tive changes within relatively short time intervals in virtu- 
ally arbitrary signals, even if the signals contain strongly pro- 
nounced periodic components, as was the case for C1-CK44. 
It should be noted that the nonstationarity factor should be 
analyzed not only for the original signal, but also for its 
smoothed (low-frequency) and fluctuation (high-frequency) 
components. 

In order to validate the proposed phenomenological sys- 
tems approach, comprehensive monitoring of seismically ac- 
tive regions such as the Kamchatka peninsula should be per- 
formed and the data should be analyzed with the FNS nonsta- 
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tionarity factor The measured characteristics should include 
geoacoustic, hydrological/hydrochemical, gasgeochemical, 
geodetic, and electromagnetic signals summarized in the in- 
troduction. 

Acknowledgements. This study was supported in part by the Rus- 
sian Foundation for Basic Research, project nos. 09-05-98543 and 
10-02-01346. 

References 

Bella, R, Biagi, R E, Caputo, M., Cozzi, E., Monica, G. D., Ermini, 
A., Gordeez, E. I., Khatkevich, Y. M., Martinelli, G., Plastino, 
W., Scandone, R., Sgrigna, V., and Zilpimiani, D.: Hydrogeo- 
chemical anomalies in Kamchatka (Russia), Phys. Chem. Earth, 
23, 921-925, 1998. 

Belyakov, A. S.: Magnetoelastic acoustic geophones for geophys- 
ical research and earthquake prediction (in Russian), Seismich- 
eskie Pribory, 33, 27-45, 2000. 

Biagi, P. E, Ermini, A., Cozzi, E., Khatkevich, Y. M., and Gordeev, 
E. I.: Hydrogeochemical precursors in Kamchatka (Russia) re- 
lated to the strongest earthquakes in 1988-1997, Nat. Hazards, 
21, 263-276, 2000. 

Biagi, P. E, Castellana, L., Minafra, A., Maggipinto, G., Maggip- 
into, T, Ermini, A., Molchanov, O., Khatkevich, Y. M., and 
Gordeev, E. I.: Groundwater chemical anomalies connected with 
the Kamchatka earthquake (M=7.1) on March 1992, Nat. Haz- 
ards Earth Syst. Sci., 6, 853-859, 2006. 

Cicerone, R. D., Ebel, J. E., and Britton, J.: A systematic compi- 
lation of earthquake precursors, Tectonophysics, 476, 371-396, 
2009. 

Descherevsky, A. V., Lukk, A. A., Sidorin, A. Y, Vstovsky, G. V., 
and Timashev, S. E: Elicker-noise spectroscopy in earthquake 
prediction research, Nat. Hazards Earth Syst. Sci., 3, 159-164, 
2003. 

Du, J., Amita, K., Ohsawa, S., Zhang, Y, Kang, C, and Ya- 
mada, M.: Experimental evidence on formation of imminent 
and short-term hydrochemical precursors for earthquakes, Appl. 
Geochem., 25, 586-592, 2010. 

Eedotov, S. A., Gusev, A. A., Shumilina, L. S., and Chemyshova, 
V. G.: The seismofocal zone of Kamchatka (in Russian), Volka- 
nologiyai Seismologiya, pp. 91-107, 1985. 

Gavrilov, V., Bogomolov, L., Morozova, Y, and A., S.: Variations 
in geoacoustic emissions in a deep borehole and its correlation 
with seismicity, Ann. Geophys.-Italy, 51, 737-753, 2008. 

Geller, R. J.: Earthquake prediction: a critical review, Geophys. J. 
Int., 131,425-450, 1997. 

Geller, R. J., Jackson, D. D., Kagan, Y. Y, and Mulargia, E: Earth- 
quakes cannot be predicted. Science, 275, 1616-1617, 1997. 

Gordienko, V., Gordienko, T, Krasnopistsev, N., Kuptsov, A., Lari- 
onov, I., Marapulets, Y, Rutenko, A., and Shevtsov, B.: Anomaly 
in high-frequency geoacoustic emission as a close earthquake 
precursor, Acoust. Phys., 54, 82-93, 2008. 

Hartmann, J. and Levy, J. K.: Hydrogeological and gasgeochemical 
earthquake precursors - A review for application, Nat. Hazards, 
34, 279-304, 2005. 

Hayakawa, M. and Timashev, S. E: An attempt to find precursors 
in the ULF geomagnetic data by means of flicker noise spec- 
troscopy. Nonlinear Proc. Geophys., 13, 255-263, 2006. 



ICEECP: Operational Earthquake Forecasting: State of Knowledge 
and Guidelines for Utilization, http://www.iaspei.0rg/d0wnl0ads/j 
Ex_Sum_v5_THJ9_A4format.pdf access: 1 October 2010, 2009. 

Ida, Y, Hayakawa, M., and Timashev, S.: Application of differ- 
ent signal analysis methods to the ULF data for the 1993 Guam 
earthquake, Nat. Hazards Earth Syst. Sci., 7, 479-484, 2007. 

Ihmle, P. and Jordan, T: Teleseismic search for slow precursors to 
large earthquakes. Science, 266, 1547-1551, 1994. 

Khatkevich, Y. and Ryabinin, G.: Geochemical an ground-water 
studies in Kamchatka in the search for earthquakes precursors 
(in Russian), Vulkanologiya i Seysmologiya, (4), 34-42, 2006. 

Kirschvink, J. L.: Earthquake prediction by animals: Evolution and 
sensory perception, B. Seismol. Soc. Am., 90, 312-323, 2000. 

Kopylova, G., Sugrobov, V., and Khatkevich, Y: Variations in 
the regime of springs and hydrogeological boreholes in the 
Petropavlovsk polygon (Kamchatka) related to earthquakes (in 
Russian), Vulkanologiya i Seysmologiya, (2), 53-70, 1994. 

Li, Y, Liu, Y, Jiang, Z., Guan, J., Yi, G., Cheng, S., Yang, B., 
Eu, T, and Wang, Z.: Behavioral change related to Wenchuan 
devastating earthquake in mice, Bioelectromagnetics, 30, 613- 
620, 2009. 

Reasenberg, P. A.: Foreshock occurrence before large earthquakes, 
J. Geophys. Res., 104, 4755-4768, 1999. 

Telesca, L. and Lovallo, M.: Non-uniform scaling features in cen- 
tral Italy seismicity: A non-linear approach in investigating seis- 
mic patterns and detection of possible earthquake precursors, 
Geophys. Res. Lett., 36, LOl 308, 2009. 

Telesca, L., Lapenna, V, Timashev, S., Vstovsky, G., and Mar- 
tinelli, G.: Flicker-noise spectroscopy: a new approach to in- 
vestigate the time dynamics of geoelectrical signals measured in 
seismic areas, Phys. Chem. Earth., 29, 389-395, 2004. 

Telesca, L., Lapenna, V, Macchiato, M., and Hattori, K.: Investigat- 
ing non-uniform scaling behavior in Ultra Low Frequency (ULF) 
earthquake-related geomagnetic signals. Earth Planet. Sci. Lett., 
268, 219-224, 2008. 

Telesca, L., Lovallo, M., Ramirez-Rojas, A., and Angulo-Brown, 
E: A nonlinear strategy to reveal seismic precursory signatures in 
earthquake-related self-potential signals, Physica A, 388, 2036- 
2040, 2009a. 

Telesca, L., Lovallo, M., Ramirez-Rojas, A., and Angulo-Brown, 
E: Scaling instability in self-potential earthquake-related signals, 
Physica A, 388, 1181 - 1186, 2009b. 

Telesca, L., Lovallo, M., and Carniel, R.: Time-dependent Fisher 
Information Measure of volcanic tremor before the 5 April 2003 
paroxysm at Stromboli volcano, Italy, J. Volcanol. and Geoth. 
Res., 195, 78-82, 2010. 

Timashev, S., Polyakov, Y, Lakeev, S., Misurkin, P., and Danilov, 
A.: Fundamentals of fluctuation metrology, Russ. J. Phys. Chem. 
A, 84, 1807-1825, 2010a. 

Timashev, S. E: Flicker noise spectroscopy and its application: In- 
formation hidden in chaotic signals, Russ. J. Electrochem., 42, 
424, 2006. 

Timashev, S. E: Fliker-Shumovaya Spektroskopiya: Informatsiya 
V khaoticheskikh signalakh (Flicker-Noise Spectroscopy: Infor- 
mation in Chaotic Signals), Fizmatlit, Moscow, 2007. 

Timashev, S. F. and Polyakov, Y. S.: Review of flicker noise spec- 
troscopy in electrochemistry, Fluct. Noise Lett., 7, R15, 2007. 

Timashev, S. F. and Vstovskii, G. V.: Flicker-noise spectroscopy 
for analyzing chaotic time series of dynamic variables: Problem 



6 



G. V. Ryabinin et al.: Earthquake precursors in the data for the Kamchatka peninsula 



of signal-to-noise relation, Russ. J. Electrochem., 39, 141-153, 

2003. 

Timashev, S. R, Polyakov, Y. S., Misurkin, P. 1., and Lakeev, S. G.: 
Anomalous diffusion as a stochastic component in the dynamics 
of complex processes, Phys. Rev. E, 81, 041 128, 2010b. 

Uyeda, S., Nagao, T., and Kamogawa, M.: Short-term earthquake 
prediction: Current status of seismo-electromagnetics, Tectono- 
physics, 470, 205-213, 2009. 

Vstovsky, G. V., Descherevsky, A. V., Lukk, A. A., Sidorin, A. Y., 
and Timashev, S. P.: Search for electric earthquake precursors by 
the method of Ricker-noise spectroscopy, Izv.-Phys. Soild Earth, 
41,513-524, 2005. 

Wyss, M. and Booth, D. C.: The lASPEl procedure for the evalu- 
ation of earthquake precursors, Geophys. J. Int., 131, 423^24, 
1997. 

Wyss, M.; Aceves, R. L., Park, S. K.; Geller, R. J., Jackson, D. D., 
Kagan, Y. Y, and Mulargia, P.: Cannot earthquakes be pre- 
dicted?. Science, 278, 487^90, 1997. 

Yokoi, S., Ikeya, M., Yagi, T, and Nagai, K.: Mouse circadian 
rhythm before the Kobe earthquake in 1995, Bioelectromagnet- 
ics, 24, 289-291, 2003. 




Fig. 1. Schematic of the measurement area (small rectangular frame 
on the left) and epicenters of largest earthquakes {Mi > 6, < 50 
km, Z) < 350 km) from 1985 to 2009, where Mi - local earth- 
quake magnitude, H - depth, D - distance from the epicenter. The 
large frame on the right shows a zoomed-in view of the positions of 
hydrogeological stations: 1 - Pinachevo, 2 - Moroznaya, 3 - Khle- 
bozavod, 4 - Verknyaya Paratunka. The solid circles denote the 
earthquakes reportedly preceded by hydrogeochemical anomalies. 
The dashed line is the axis of the deep-sea trench. The earthquakes 
were selected using the catalog of Geophysical Survey, Kamchatka 
Branch of the Russian Academy of Sciences. 



G. V. Ryabinin et al.: Earthquake precursors in the data for the Kamchatka peninsula 



7 




Q 300 



91 92 93 94 95 96 97 98 99 00 01 02 03 04 05 OS 07 06 09 



f (years) 



Fig. 2. Comparison of nonstationarity factor Cj (T=600 days, AT=3 days) for the GK-1 chlorine-ion concentration time series with seismic 
activity: Vq - source signal; Co - nonstationarity factor for Vg, Cr - nonstationarity factor for the low-frequency component of Vg, Cf 
- nonstationarity factor for the high-frequency component of Vg, Mi - local earthquake magnitude, D - distance from the epicenter. Solid 
triangles denote sample Cr spikes preceding large earthquakes. Crosses denote sample Cr spikes not related to large seismic events. 




Fig. 3. Comparison of nonstationarity factor Cj (T=600 days, Ar=3 days) for the GK-44 chlorine-ion concentration time series with 
seismic activity: Nomenclature as in Fig. 2. Solid triangles denote sample Cp spikes preceding large earthquakes. Crosses denote sample 
Cf spikes not related to large seismic events. 
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Fig. 4. Nonstationarity factors for GK-1 and GK-44 chlorine-ion concentrations and Z160 G-1 geoacoustic emissions in the time frame 
around the 8/10/2001 earthquake. Mi - local earthquake magnitude, D - distance from the epicenter. The double-headed arrows denote the 
time intervals between the nonstationarity factor spikes and earthquake itself. 



